The gluon propagator close to criticality 
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There are good reasons that the deconfinement phase transition of pure Yang-Mills theory at finite 
temperature should also be reflected in the behavior of gauge-fixed gluonic correlation functions. 
Understanding this in detail would provide another important example of how physical informa- 
tion can be extracted from gauge- dependent correlations, which is not always obvious. Therefore, 
herein we study the behavior of the Landau-gauge gluon propagator of pure SU{2) across the phase 
transition in 2-1-1 and 3-1-1 dimensions in order to assess to what extend the corresponding criti- 
cal behavior is reflected in these correlations. We discuss why it should emerge from a continuum 
perspective and test our expectations in lattice simulations. A comparison with SU{3) furthermore 
reveals quite clear indications for a sensitivity of the gluon propagator to the order of the transition. 

PACS numbers: 12.38.Aw,12.38.Lg,11.15.Ha,12.38.Mh,25.75.Nq 
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I. INTRODUCTION 

Our knowledge of even the most basic characteristic 
features of the QCD phase diagram is still rather uncer- 
tain. A rich structure of possible phases and transitions 
between them is being discussed in various regions, see, 
e.g. [H-Q- Understanding how these phase transitions 
manifest themselves both theoretically and experimen- 
tally in a joint effort is the great challenge of strongly 
interacting matter research. 

One of the well established properties of QCD with 
2-1-1 flavors of dynamical quarks is a crossover at vanish- 
ing chemical potential in a temperature range of about 
150-180 MeV, e.g. i^. However, many observables 
nevertheless appear to reflect the deconfinement phase 
transition that would occur if all quarks were sufficiently 
heavy. In fact, whether and where in the so-called 
Columbia plot the change from a genuine phase transi- 
tion to a crossover occurs might still change when going 
beyond QCD as an isolated theory. In particular, if one 
also considers the quarks' fractional electric charges and 
couplings to more than one gauge group as in the Stan- 
dard Model, a center-like symmetry reemerges which can 
break at the transition as in the quenched case Q. 

This finite-temperature deconfinement phase transi- 
tion of the pure SU{N) Yang-Mills theory in d -f 1 di- 
mensions follows the Zn symmetry breaking pattern of 
a d-dimensional g-state Potts model with q = N. It is 
of second order for N = 2 (Ising universality class) in 
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d = 2 and 3 and for iV = 3 in d = 2. The self-duality of 
the 2-dimensional Potts models thereby manifests itself 
in the free energies of the confining electric fluxes being 
mirror images around criticality of those of center-vortex 
ensembles with twisted boundary conditions [9-11]. 

Another tool to determine the dynamics of order pa- 
rameters for the Yang-Mills phase transition, in particu- 
lar the Polyakov loop, the chiral condensate and to some 
extent thermodynamic bulk properties, have been gauge- 
fixed correlation functions, in particular in the Landau 
gauge. While it has been known since decades how 
the chiral phase transition and in particular its criti- 
cal properties manifest themselves directly in the gauge- 
fixed quark correlation functions, it has only recently be- 
come clear how to compute the non-perturbative effective 
potential of deconfinement order parameters from the 
gauge-fixed propagators 7, 12-18]. These analyses entail 
that already the gauge-fixed propagators carry the crit- 
ical properties of the confinement-deconfinement phase 
transition, which has been used explicitly in [13|. These 
findings have also been studied using lattice simulations, 
see in particular [3 113 j though systematic errors are 
still a matter of ongoing research [18l422j . 

In the present paper we extend these investigations. 
We first review in Section [ll] how critical properties of 
the phase transition are incorporated in the gauge-fixed 
correlation functions. We argue that the gluon propaga- 
tor should be directly sensitive to the phase transition 
and, in particular, that it should refiect critical scaling. 
In Section IIIII the behavior expected from these argu- 
ments is compared to results from lattice simulations in 
three and four dimensions for SU{2) Yang-Mills theory. 
In both cases the phase transitions are of second order, 
but one should be able to distinguish the different critical 
exponents corresponding the 2d and 3d Ising universal- 
ity classes, respectively. We also briefly revisit the flrst- 
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order case of four-dimensional SU{3) Yang-Mills theory 
to assess the sensitivity of the gluon propagator to the 
order of the transition in Section [IVI Our summary and 
conclusion are given in Section [V] Some additional re- 
sults on the momentum dependence of the propagators 
are deferred to Appendix I3 while an extended discus- 
sion of the systematic and statistical errors is given in 
Appendices |B] and O respectively. 



Usually /? is bound by its mean field value, i.e. /3 < 1/2, 
and the derivative of M{t) with respect to t therefore 
diverges at the critical temperature. 

A natural order parameter for confinement is provided 
by the Polyakov loop. 



L(x) = -5-trP exp 
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l/T 



dr Ao{t, x) 



(5) 



II. CRITICAL BEHAVIOR AND THE GLUON 
PROPAGATOR 

The second-order finite-temperature deconfinement 
phase transitions of SU{2) Yang-Mills theory in four, and 
both SU{2) and SU{3) in three dimensions, are charac- 
terized by correlation lengths which diverge at the criti- 
cal temperature in the infinite volume limit. As a conse- 
quence, only the long-range properties of the theory mat- 
ter near criticality. These are determined by symmetries 
and dimensionality of the system, and the concepts of 
universality and scaling apply. The second-order SU {N) 
cases belong to Potts/Ising universality classes in two 
and three dimensions as mentioned above. In particular, 
the SU{2) Yang-Mills theories in three and four dimen- 
sions that we are mainly interested in here, belong to 
the classes of the two-dimensional and three-dimensional 
Ising models, respectively. 

Critical behavior manifests itself in singularities of 
thermodynamic functions at the critical temperature Tc 
(and vanishing external field). In particular, these singu- 
larities show as characteristic non-integer powers of the 
reduced temperature 
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(1) 



The critical exponents in these power laws only depend 
on symmetry-breaking pattern and dimension, and are 
thus a characteristic feature of a universality class. Scal- 
ing and hyper-scaling (in four dimensions or less) entail 
that all critical exponents can be expressed in terms of 
two independent ones, say v and rj. The first is the crit- 
ical exponent of the diverging correlation length ^, i.e.. 



(2) 



close to a second-order phase transition. The second 
is defined from the (in general renormalization-group- 
dependent) connected propagator G{p) of the order pa- 
rameter, whose zero momentum contribution defines the 
susceptibility which diverges with exponent 7, 



G(0) cx l^r'' , with 7 = i^(2-7]). 



(3) 



In the broken phase, an order parameter AI scales with 
the exponent /?, 
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It depends on the zero component Aq of the gauge field, 
and its correlation length is inversely related to the string 
tension, i.e., = T/a{T) for T < T^. Therefore, with 
second-order transition, a cx {—ty for t — > 0^.^ 

In a gauge-fixed approach the theory is described in 
terms of, in general, gauge-dependent correlation func- 
tions, i.e. A^-point functions of the gauge field. These cor- 
relation functions have to contain the information about 
the physics at criticality, even though they change under 
gauge transformations. 

Indeed it has been shown how to extract the expecta- 
tion value of the (logarithm of the) Polyakov loop from 
the gauge-fixed propagators in general gauges in terms of 
the full effective potential Kff of {Aq) |13 - [14| . The value 
of {Aq) at the minimum of Vcs directly relates to the ex- 
pectation values of the eigenvalues of the (logarithm of 
the) Polyakov loop, it is gauge invariant and also provides 
an order parameter. Other recent examples are dual con- 
densates [23 - [28j , which can likewise be computed from 
gauge-fixed correlation functions [isl - fTsl . [2§| , and in the 
quenched case determine the Polyakov loop directly. 

In [1^ [iJl the potential has been computed in Landau 
gauge; the critical temperature agrees well with lattice 
results. In [l3| the same result was obtained within a 
computation in Polyakov gauge, a necessary requirement 
for the results to be gauge-independent as they must be. 
Moreover it has been verified that in SU{2) the poten- 
tial and hence the order parameter does indeed show the 
scaling of the Ising universality class. 

An interesting consequence of these computations and 
the related formal analysis is that the gauge-fixed prop- 
agators necessarily also show critical scalin g. Here we 
briefiy review the corresponding results in |12l4l4| . In 
(1^ it has been shown that the effective potential of {Aq) 
in a general gauge can be computed solely from the ghost 
and gluon propagators in the constant background Aq: 



Tr {C^C').dkR'C , (6) 



M{t)(x\tf with (i = -v{d - 2 + -q) 



(4) 



^ Above Tc, the Polyakov- loop correlation length ^+ is related 
to the dual string or interface tension tr oc (5+)^"'*, a well- 
defined order parameter for center-symmetry breaking [23j . With 
second-order transition, a oc t^^ for t — > 0+, and (hyper-scaling) 
fj, = {d — 1)1/. 
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where k is an infrared cut-off scale implemented by the 
regulator i?^'' in the functional renormalization group 
equations. At large cut-off scales the effective potential 
vanishes. Hence, integrating the above equation ^ from 
the trivial potential at A; — >■ oo to fc = provides us with 
Vcs{Ao) = K!ff,k=o(^o) the full non-pcrturbative effec- 
tive potential of (^o)- 

Since Vcs{A-q) is the potential of the order parameter 
(log L) , it has to show critical scaling (if the transition 
is of second order). Consequently, also the propagators 
will be sensitive to critical scaling, in general. 

The above setting even allows us to extract the scaling 
of the longitudinal gluon. Its propagator schematically 
is of the form, 
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{~DI+jF)Zl 



(7) 



where 11^ stands for the projection operator longitu- 
dinal to the heatbath in the presence of the back- 
ground ^0, and iDq = po + gA^- _The prefactor Zl = 
Zl{—Dq,jF\Ao) depends on ^o, and separately 
and reduces to the well-known longitudinal (inverse) Lan- 
dau gauge dressing function for Aq — 0. 

The solution of the equations of motion for Aq is 
the order parameter related to (log L) and derives from 



9Ao 



0. 



(8) 



Ao=ao 



Hence it scales according to (HJ , and the effective poten- 
tial also shows critical scaling, i.e. 



dAl 



(ao)(x lip. 



(9) 



with 7 = 1^(2— 77), see ([3]). Up to possible dependencies on 
the background field, the inverse longitudinal gluon prop- 
agator is proportional to {AlAl)~^ oc d'^Vcfi/ dA-^. 
In other words, the mass gap is given by the curvature of 
Kff(^o) in (El), and hence determines the singular behav- 
ior of the gluon propagator at vanishing frequency and 
momentum, i.e. 



(10) 



We conclude that at least the longitudinal or electric 
propagator ([7]) scales as the propagator of an order pa- 
rameter, see (O. In [13| this has in fact been used to 
determine the critical scaling directly from the longitu- 
dinal gluon propagator. 

Note however, that a similar analysis cannot be made 
for transverse gluons and ghosts. This is due to the fact 
that their inverse propagators contain no ylg-derivatives 
and hence no dependence on the curvature term in 
Consequently, we cannot infer from this simple argument 
whether one should see critical scaling in the transverse 
gluon and ghost propagators or not. However, it is sug- 
gestive that it at best arises in sub-leading order, that is. 



that singularities might appear in higher T-derivatives. 
This is supported by lattice results, see [T8l - [20l Isoj and 
appendix [Cj 

Now we turn to lattice Landau gauge and, in partic- 
ular, discuss the differences of simulations there to the 
background Landau gauge scenario used above. The 
screening mass of a two-point correlation function, the 
propagator D{p), defined as 



(11) 



is a candidate for critical behavior. Note that the corre- 
sponding pole mass [30], if it exists, will in general only 
be sensitive to correlation length s up to its inverse size. 
Since all results available so far [SOj seem to exclude a 
zero pole mass of the gluons, it is therefore to be expected 
that the pole mass will not exhibit critical behavior. 

The temperature-dependence of the electric and mag- 
netic gluon propagators and the ghost propagator in lat- 
tice Landau gauge has been investigated previously, for 
details see [l8l - l22i l3ll - [39| . One observes that only the 
electric propagator is sensitive to the phase transition 
[iM[3[23- However, the available lattice data do not yet 
allow to draw firm conclusions about a possible scaling 
behavior at criticality. Nevertheless, the susceptibility 
defined as the temperature derivative of its screening 
mass M^, 



X 



dT 



is sufficiently sensitive to determine the phase transition 
temperature within the systematic errors, see [l^ and 
Appendix |B] 

The aim of the present investigation is the further as- 
sessment of the critical behavior. This will be performed 
using lattice simulations in three and four dimensions 
for the gauge group SU(2). In four dimensions, for this 
purpose, a significantly refined temperature mesh and 
increased statistics were used compared to the previous 
study ■ The three-dimensional system, investigated 
here for the first time, has the advantage, compared to 
the four-dimensional one, that the expected critical ex- 
ponents are larger, and therefore the signal might be 
stronger. 

In lattice Landau gauge, naively, one would expect a 
similar behavior as for the background Landau gauge dis- 
cussed above. However, the two differ in at least one im- 
portant qualitative aspect. All currently employed prac- 
tical implementations of the lattice Landau gauge, in- 
cluding the one used herein, exhibit a screening mass. 
Given its non-perturbative origin and appearance in a 
gauge-fixed correlator, it will be referred to as a Gribov 
mass [40,]^, m,gi-ibov This mass is present at all temper- 



^ Note that in the original work [40| | the so defined screening mass 
was infinite. 
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Model of the electric screening mass 
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FIG. 1: Cartoon of the electric screening mass and suscep- 
tibilities of the form (|12|) and (|13|l . with critical exponents 
7/2 (divergent case) and 7 (finite case) with 7 = 1.2372(5), i. 
e., with the 3-dimensional Ising universality class. The non- 
universal constants are mgribov ~ 0.2 GeV, a = —1/4 GeV 
and 6 = 3/4 GeV. 



atures and hence interferes with critical scaling. This 
entails that the inverse propagator can still have a mass- 
gap at criticality. It is important to emphasize that the 
effective potential V^s for the order parameter ao in such 
a gauge still satisfies (|6]). Hence, the propagators must 
reflect critical scaling nevertheless. This implies that the 
screening mass Mf of the electric propagator in ([TT]) con- 
tains both, a regular Gribov mass contribution as well as 
a singular critical one. From Eq. ([6]) one concludes that 
the critical contribution is determined by the exponent 7. 
Unfortunately, this general argument is not sufficient to 
disentangle the two contributions in a unique way. The 



two perhaps most natural candidates to explore at least 
as approximate descriptions near criticality might be to 
either simply add the two mass contributions directly, i.e. 

Mf (t) = mgribov + a± \tyl^ , (12) 

where a-|_ = a is a non-universal coefficient for temper- 
ature i > and = b for temperatures t < 0; or, 
alternatively, add two corresponding self-energy contri- 
butions to the inverse propagator which would lead to a 
sum of squares for the total mass, i.e. 

Mf (t) = y^m^^ibov + 4 1^1^ ■ (13) 

A simply testable criterion to at least discriminate these 
two cases is provided by the corresponding susceptibil- 
ities which near criticality in the first case (IT^ be- 
haves as lil'^/^'S while with (US]), x^ 1*^"^ For 
1 < 7 < 2, in the 2d and 3d Ising universality classes 
one has 7 — 7/4 — 1.75 and 7 = 1.2372(5), respec- 
tively [4l| . This implies that (|12p leads to a susceptibility 
X^ = dtMf which would diverge at the critical tempera- 
ture (in the infinite volume limit), while (jl3p would lead 
to a vanishing one. This is illustrated by a cartoon of the 
four-dimensional case in Figure [TJ 

While there will of course neither be a divergence nor 
a strict zero in x^ in ^ finite volume, the lattice results 
clearly favor an enhancement of x^ near Tc correspond- 
ing to (jl2p rather than a suppression, as can be seen from 
comparing the cartoon in Figure [1] with Figures [2] and |4] 
below. So we conclude that ([TS]) can relatively safely be 
ruled out. A more conclusive study of the expected be- 
havior based on a self-consistency analysis of ([6|) would 
be clearly desirable but will be postponed to future work. 

Furthermore, because the Gribov mass itself is signif- 
icantly volume dependent [sO, 42- 45], it is not possible 
at present to clearly identify the critical behavior in the 
electric screening mass or its temperature derivative from 
a finite-size-scaling analysis of lattice results. The Gri- 
bov mass creates a mass gap whose volume dependence 
especially at criticality is yet to be determined. A signifi- 
cant amount of different volumes, and possibly discretiza- 
tions, will be necessary to disentangle the two intertwined 
finite- volume effects on the electric screening mass, from 
critical scaling and from the Gribov mass, requiring sub- 
stantially more resources than presently available to us. 

Finally, it should be noted that the ansiitze ([T2|) and 
(HH) are expected to be only valid in the critical region. 
The Gribov mass in these fits will therefore not coincide 
with the one at zero temperature, as further, sub-leading 
temperature-dependent contributions are neglected. 



III. LATTICE RESULTS 

The gluon propagator is determined using standard 
lattice methods, for details see [l^E^. We use the min- 
imal Landau gauge [3Ci| with the implications discussed 



in Sectionini In particular, the gluon propagator exhibits 
a finite Gribov screening mass already at zero tempera- 
ture, while the ghost behaves essentially tree-level-like 
[SmI, I47I, [48| . In four dimensions, in addition to the 
data from [l8| , further results using a finer discretization 
at the price of smaller physical volumes, are used. In 
three dimensions, two different discretizations are used 
as well. For the main purpose of this work only the elec- 
tric screening mass is of importance. For completeness, 
also the finite momentum results are reported, but rele- 
gated to Appendix El 

For setting the scale a string tension of (440 MeV)^ has 
been used, based on the a values determined in [4§| 
for three dimensions and in [50| for four dimensions. This 
particular value was chosen to ensure comparability to 
the previous investigations [3, [l^ . The critical /3 values 
are taken from [51] for three dimensions and from [SQ] for 
four dimensions. A discussion of systematic effects due 
to the scale setting can be found in Appendix IB 11 



A. Three dimensions 

The results for the electric screening mass, as well as 
its susceptibility in three dimensions are shown in Figure 
121 It is clearly visible that the maximum of the suscepti- 
bility marks the critical temperature, which has been de- 
termined independently using analyses of the free energy 
and quantities derived from it |50l [sH . Furthermore, a 
sizable volume scaling is seen for the smaller spatial vol- 
umes with Nf ~ A time slices. For the largest spatial 
extent the peak in the susceptibility only becomes some- 
what sharper, but not significantly higher. This might 
be due to the mixing with the Gribov mass as discussed 
in Section |lTl Note, however, that the singular contri- 
bution to = dtMf (X from our model 
would predict a finite-size L scaling of the maximum as 
xSax oc L(i-7/2)/>' = ii/8 for three-dimensional 5t/(2), 
corresponding to a factor 2^/^ or about only 9% increase 
with doubling the spatial lattice size which is still within 
errors. Furthermore, increasing Nt from 4 to 6, which 
effectively reduces the physical L for equal spatial lattice 
sizes by 2/3, we observe that the results remain qualita- 
tively unchanged, at least for the two values checked here. 
The details might still be influenced by lattice artifacts. 

We reiterate that the screening mass, which is 
(semi)positive by definition, is non-zero at the phase 
transition. Moreover, because the peak in its tempera- 
ture derivative marks the phase transition, the screen- 
ing mass does not have a minimum there, but a point of 
inflection.^ It is also well established that there is a finite 



Electric screening mass in three dimensions! 



^ From the coarser temperature resolution in IT^ it was suggested 
that the transition occurred at the maximum of the electric gluon 
propagator at zero momentum, i.e. for minimal Mf, which is 
slightly off. However, this merely implies that the phase transi- 
tion temperature should be set to the peak in the susceptibility 
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FIG. 2: The electric screening mass (top panel and zoomed in 
in the phase transition region in the middle panel) as a func- 
tion of the reduced temperature for various lattice sizes. The 
bottom panel shows the corresponding susceptibility, where 
some temperature points have been dropped to reduce the 
statistical errors. 
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screening mass in minimal Landau gauge at zero tem- 
perature 4J|. The possibility that it remains constant 
throughout all temperatures up to the phase transition, 
and that it starts to increase abruptly there, would im- 
ply that the non-universal coefficient a_ = & = and 
is included here. For our fits we thus assume a con- 
stant contribution corresponding to a Gribov mass intro- 
duced by the gauge-fixing, which should be temperature- 
independent, to which we add a form resembling a criti- 
cal behavior as in Subleading contributions around 
criticality will be neglected, leading to the Ansatz for the 
screening mass. 



Fit of the electric screening mass | 



0.4 



M^{t) = TOgribov + 0{t)at— + 9{~t) b {-ty 



(14) 



where we have furthermore introduced independent ex- 
ponents 7+ and 7_ above and below Tc in order to assess 
to what extend the results are consistent with a unique 
7+ = 7__. As a side remark note that the Ansatz (fT4|) 
for small t in principle also includes the form in for 
which one should obtain 7+ = 7_ = 27. 

Performing the fits in the reduced-temperature interval 
[—0.1, 0.1] yields the results shown in tableHl An example 
of the quality of the fit is shown in Figure |31 There are 
a number of observations: 

First, the fit form (fT4|) describes the screening mass 
rather accurately, and within the statistical errors also 
the susceptibility in a satisfactory way. The resulting 
exponents are all in the range 1 < 7 < 2 (and thus in- 
consistent with ([T3l) ). but they do not agree with one 
another very well, which is only slightly improved if the 
fit interval is further reduced, assuming a smaller critical 
region. 

Given that the fits provide a rather good description 
of both the screening mass and the susceptibility there 
are a number of possible reasons for the observed non- 
singularity in the data itself. The first is that the volumes 
are just too small and the discretization too coarse to be 
close enough to a critical behavior. Given that the crit- 
ical behavior is well visible on comparable lattices for 
other observables [5l| , this appears at first glance to be a 
rather unlikely explanation. However, it is possible that 
the electric gluon propagator lags significantly in the de- 
velopment of the critical behavior behind those other ob- 
servables. Indeed, the gluon propagator requires already 
at zero temperature very large volumes to reach even 
its qualitative infinite- volume behavior |4 2[-[45| . This ap- 
pears to be true also at finite temperature [i3,[2l|. There- 
fore this appears to be a possible explanation. 

The second possibility is that the critical region is very 
narrow. Then it would be necessary to investigate the 
region around the critical temperature with a much finer 
mesh. In this case also significantly enhanced statistics 
would be necessary, since the numerical derivative for 
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instead. Otherwise the results of [l^ remain unchanged. 



FIG. 3: The quality of the fits compared to the data from the 
4 X 128^ lattice. The standard fit uses the values from Table 
in As an alternative, a fit using only a more restricted reduced 
temperature interval [—0.05, 0.05] is also shown (dotted line), 
which yields 7_ = 1.52'^l'^ and 7+ = 1.61^^4. Fits have been 
performed on the screening mass using all data points, and 
the derivatives have then been calculated analytically. The 
error bands are also shown. 



the determination of the susceptibility requires that the 
statistical error must be much smaller than the distance 
between two points. 

Finally, it is of course still possible that the electric 
gluon propagator simply does not show critical behav- 
ior. At the moment this cannot entirely be excluded, 
despite the general continuum arguments in favor given 
in Section |TT1 This could be due, e.g. to a more compli- 
cated interference with the Gribov mass than the two 
possibilities considered in and (IT51) . An alternative 
possibility will be discussed in Section UlI CI below. How- 
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TABLE I: The fit parameters for tlie fit (|14|) for tlie tliree-dimensional case. For tlie determination and placement of the 
statistical errors see Appendix[Cl Only the change in the last digits is indicated. See Appendix[B]for a discussion of systematic 
errors. 



Lattice 


m-gHbov [GeV] 


a [GeV] 


b [GeV] 


7+ 


7- 


4 X 32^ 


0.37019:^3 


0.85^^ 


-0.78tl 


i.88;l 


1.8i;l 


4 X 64^ 


0.2637;^ 


0.90;* 


-0.67+3 


1.68;2 


1.63^2' 


4 X 128^ 


0.221;^ 


0.735"? 


-0.368+i^ 




1.134+^ 


4 X 256^ 


0.218;3 


o.73;l 


-0.33^3 


1.39;i^ 


1.02;^ 


6 X 128^ 


0.2571^5 


o.8o;^ 


-0.54+^ 


1.63+;6^ 


1.38+1 



TABLE 11: As Table HI but for four dimensions. 



Lattice 


mgribov [GeV] 


a [GeV] 


b [GeV] 


7+ 


7- 


4 X 46^ 


0.251^ 


l-16r-6 


-o.o4ir 


1.26;iJ 




6 X 48^ 


0.25;^ 


1.5li 


-o.o7+;f 


l-54o:S^ 


o.e+f 



ever, since the behavior of the Polyakov loop can be de- 
rived from the propagators TSj, this appears somewhat 
unlikely, though by far not impossible. 



B. Four dimensions 

Essentially the same qualitative behavior, though 
within a different universality class, is expected in four 
dimensions. The corresponding electric screening mass 
and its susceptibility are shown in Figure U) The quali- 
tative picture is the same as in three dimensions. Note, 
however, that because the critical temperature is about 
one third lower in four as compared to three dimensions, 
the physical volumes are about 50% larger here than they 
would be for the same number of lattice points in three 
dimensions. Even though our volumes are nevertheless 
still smaller than in the three-dimensional case, the sus- 
ceptibility is significantly more peaked, showing that the 
phase transition leaves a much more significant imprint. 
However, the peak of the susceptibility is somewhat dis- 
placed as compared to the phase transition point. One 
possible reason could be systematic errors in the scale set- 
ting. Other, equally possible, reasons may be statistical 
fluctuations, finite volume effects or discretization effects. 
Permitting even only a systematic error of a percent , the 
correct phase transition temperature is obtained within 
the total error. Various such sources of systematic un- 
certainties are further explored in Appendix |B] 

It is possible to perform the same fits as for the three- 
dimensional case with the Ansatz The results are 
given in Table HIl and shown in Figure [S] The situation is 
qualitatively the same as in three dimensions. The result- 
ing critical exponents tend to be somewhat smaller than 
in three dimensions as expected from universality. The 
temperature dependence of the electric screening mass 
below Tc is weaker than before, and in is fact within er- 
rors compatible with being esssentially constant. With 
critical scaling, the prediction for its maximal slope at 
criticality would be roughly Xmax L^'^, which would 



Phase transition vicinity in four dimensions 




-0.05 0.05 0.1 

t 



Susceptibility | 




FIG. 4: The electric screening mass (top panel) as a function 
of the reduced temperature for two different lattice sizes. The 
bottom panel shows the corresponding susceptibility, where 
some temperature points have been dropped to reduce the 
statistical errors. 
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Fit to the electric screening mass | 



Fluctuation-supressed screening mass | 
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FIG. 5: The quality of the fits compared to the data from the 
4 X 46'^ lattice. The standard fit uses the values from Tabled 
As an alternative, a fit using only a more restricted reduced 
temperature interval [—0.05, 0.05] is also shown (dashed line), 
which yields 7_ = 0.5!j;g''^ and 7+ = 1.5^2- Fits have been 
performed on the screening mass using all data points, and 
the derivatives have then been calculated analytically. 



predict a 25% difference between the two volumes used 
here. The somewhat stronger volume dependence of the 
maximum as compared to the three dimensional case 
might be reflected in Figure S] but this is well within the 
present errors still, and thus not significant. 

A further point to consider in the four-dimensional 
case, in contrast to the three-dimensional one, is renor- 
malization. In three dimensions the gluon propagator 
is finite, and no renormalization is necessary. This is 
not the case in four dimensions. As a consequence, 
the screening mass, in contrast to a pole mass, is 
renormalization-scale dependent (30| . Here, the renor- 




FIG. 6: The fluctuation of the screening mass (|I5|) in four 
dimensions from the 4 x 46"' lattice. 



malization is performed by requiring fi^D^dijT) = 1 at 
/i = 2 GeV. However, this is problematic for two rea- 
sons. First, this introduces a temperature dependence 
of the renormalization constant. Secondly, the renormal- 
ization prescription for the gluon propagator in Landau 
gauge is linked to the one for the ghost propagator which 
might lead to a potential inconsistency, since also the 
ghost propagator can acquire a temperature-dependent 
renormalization. In the present case, however, these ef- 
fects can be estimated and are found to be well within 
the statistical errors, see Appendix IB 21 



C. Alternative observables 

Since in both three and four dimensions the gluon 
propagator does not show a behavior which can immedi- 
ately be interpreted as a pure critical behavior, a natural 
question might be whether one can suppress the influence 
of non-thermal contributions depending on whether they 
are due to discretization, finite volume artifacts, genuine 
quantum effects, or the Gribov mass as mentioned. 

For this purpose the fluctuations of the screening mass 
can be used, which are defined as 



1 



D^{0) = < Z?(0) > - (/d(0))' 



(15) 



This essentially amounts to determining the pure fluctu- 
ating part. 

The resulting screening mass for the four-dimensional 
case is shown in Figure [B] There are a number of obser- 
vations which can be made, despite the significantly in- 
creased statistical uncertainty. First, the global features 
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of the screening mass are maintained. Secondly, even af- 
ter removing the non-fluctuating part there is a non-zero 
baseline, which implies that the Gribov mass from the 
gauge fixing is present even in the fluctuations, though, 
of course its infinite volume behavior cannot be inferred 
from these results. Thirdly, the slope again appears to 
increase at the phase transition, which is not a statisti- 
cally very sound result at present, however. Nonetheless, 
the behavior is suggestive, and this would imply a similar 
phase transition signal also in the fluctuations. 

In relation to the original goal of isolating critical be- 
havior, however, there is no sign of improvement as com- 
pared to the screening mass. 

A second alternative would be the magnetic screening 
mass. From the results in [31 and the results shown in 
Appendices [X] and O it can already be inferred that it 
will not exhibit a strong signal of the phase transition. 
Indeed, a more detailed investigation of its susceptibility 
shows that, if there is any signal of the phase transition 
encoded in the magnetic screening mass at all, it is cer- 
tainly substantially weaker than that in the electric one. 
The same applies to the gauge-fixing sector in form of 
the Faddeev-Popov ghosts. 



IV. A NOTE ON FOUR-DIMENSIONAL SU{3) 
YANG-MILLS THEORY 



The physically more relevant case of four-dimensional 
SU{3) Yang-Mills theory exhibits a (weak) first-order 
transition [521 ■ It h as p reviously been studied with func- 
tional methods [iMl IMM M, and on the lat- 
tice [l^ [l^, [S^]. The correlation length remains finite, 
and one does not expect critical scaling. In principle, 
one should be able to study the finite-size scaling with 
integer exponents characteristic of first-order transitions 
which is rather expensive to investigate in lattice simu- 
lations, however. 

Here we analyze the data from H^l supplemented 
by additionally generated temperature points and with 
improved statistics in a closer window around the tran- 
sition for comparison. This leads to an electric screening 
mass and a susceptibility in four-dimensional SU{5) as 
shown in Figure [T] There is no significant temperature 
dependence below Tc in the range considered here. The 
constant value serves as the baseline in our fit. Above 
Tc we allow a temperature dependence of a form which 
is typical for a first-order transition in addition to the 
constant contribution. With these assumptions we ob- 
serve that the data in the reduced temperature interval 
[—0.1,0.1] is well described by a fit of the form 



Phase transition vicinity | 



Aff = mgribov + 0(t) aVsTt, 



(16) 



with parameters mgi-ibov = 0.222^5 GeV, a — 0-7^f GeV, 
and S = 0.06^37. 

The data provides rather clear indication that a pro- 
nounced jump forms at the transition characteristic of its 
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FIG. 7: The electric screening mass and susceptibility for 
four-dimensional 5*17(3) Yang-Mills theory on 4 x 34^ lattices 
close to the phase transition, compared to a fit of type (|16[) . 



first-order nature. A much more detailed analysis espe- 
cially of the finite volume behavior of these results will 
be necessary to substantiate whether a discontinuity in 
the electric screening mass develops in the infinite volume 
limit or not. The fact that the transition is only weakly 
first order makes such an analysis particularly demand- 
ing. At present we can conclude, however, that the be- 
havior looks significantly different from the second-order 
cases discussed in Section IIIII An emerging discontinu- 
ity at the phase transition appears to be indicative of 
a sensitivity of the longitudinal gluon propagator to the 
order of the transition. One should nevertheless be wary 
of lattice artifacts ^] in this exploratory study also. 
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V. CONCLUSIONS AND OUTLOOK 

In summary, we argue that critical behavior of 
the confinement-deconfinement phase transition of pure 
Yang-Mills theory should be reflected in the longitudinal 
gluon propagator, in particular, in the electric screen- 
ing mass. Our data shows in all cases, for SU{2) and 
SU{5), that the phase transition can be clearly identi- 
fied by a peak in its temperature derivative. Although 
it turns out to be rather difficult to disentangle a sin- 
gular critical contribution to the electric screening mass 
from a regular one, which is presumably predominantly a 
Gribov copy effect, it seems that there is such a contribu- 
tion in SU{2), and that this contribution changes when 
going from three to four dimensions in a way consistent 
with the change from the 2d to the 3d Ising universal- 
ity class. Our results for the critical exponents are still 
somewhat off, however, which might be mainly due to 
an insufficient way of isolating the corresponding con- 
tributions. Whether our qualitative findings will finally 
be backed by a more quantitative analysis on the ba- 
sis of universality arguments, including a verification of 
the expected finite-size scaling, as it is well established 
for various other observables, will have to await further 
study. 

In comparison, the electric screening mass of pure 
SU{3) in four dimensions appears to show the disconti- 
nuity indicative of the first-order phase transition there. 
This seems to indicate that the electric gluon propagator 
is sensitive to the order of the transition, likewise. 

To ultimately clarify all this will require both substan- 
tially larger spatial volumes as well as finer discretiza- 
tions. It will be important to further pursue this for two 
reasons: One is that identifying critical behavior in this 
simplest of observables will improve our understanding of 
the driving mechanisms behind critical behavior and the 
order of the phase transition in general. Secondly, this 
clarification is of great interest in the construction of re- 
alistic truncation schemes for functional methods which 
have the potential of providing a genuine first-principle 
approach to describe full QCD with various quark flavors 
and masses at finite temperature and density • 

Thus, to further pursue these investigations is of prime 
importance in our efforts to combine both numerical 
lattice and functional continuum methods for a reliable 
determination of the QCD phase diagram. 
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Appendix A: Full momentum dependence 

Besides the critical behavior manifest at zero momen- 
tum also the full momentum dependence of the propa- 
gators is important, both for combining with functional 
methods as well as to estimate the reliability of other 
approximation techniques, like e. g. hard-thermal-loop 
calculations close to the phase transition. 

Hence, in this appendix the full momentum depen- 
dence of both parts of the propagator, the electric and the 
magnetic one, for the zeroth Matsubara frequency will be 
presented. Higher Matsubara frequencies can then be ap- 
proximated rather accurately with this information jl8l | . 
Results for the Faddeev-Popov ghost can be found else- 
where [ilIilllO,!!!. 

The results in three dimensions are shown in figure [5] 
for the 4 x 128^ lattice. For four dimensions the results 
for the 6 x 48^ lattice are shown in figure [9l The re- 
sults for the 4 x 46'^ lattice can be found in [l^. For 
both dimensions, the results are qualitatively similar, up 
to the differences present already at zero temperature. 
In particular, the abrupt change in the electric screen- 
ing mass discussed in the main text finds its echo at 
finite momentum. The magnetic propagator shows no 
pronounced dependence on the temperature, except for 
a general suppression. 

It is an interesting observation that the smoothness of 
the electric screening mass in three dimensions compared 
to four dimensions also surfaces in the finite momentum 
behavior. This again can be interpreted as a sign that 
Yang-Mills theory in four dimensions is 'closer' to a first 
order transition than three dimensions, as is known from 
other observables. 
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Appendix B: Systematic errors 

Calculations like the present ones are sensitive to var- 
ious systematic errors, as has been pointed out repeat- 
edly 18-22]. In the following various such errors will be 
discussed in turn, including the scale setting, renormal- 
ization, finite-volume effects and aspect ratio issues. 
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Magnetic gluon propagator in three dimensions 




Magnetic gluon propagator In four dimensions 




FIG. 8: The full momentum dependence of the zero-mode 
electric (top panel) and magnetic (bottom panel) gluon prop- 
agator in three dimensions on a 4 x 128'^ lattice. 



1. Scale setting 

The scale in the present calculations has been set in 
four dimensions using an interpolation formula of the 
type, 

r- Ci C2 C3 C4 
"^=^ + ^ + ^ + ^- 

with cr, the string tension, set to (440MeV)^, and where 
the coefficients have been fitted to the results from [50l |. 
yielding ci = 31.9, ci = -237, C3 = 574, and C4 = -444. 
However, since reduced temperatures of size 0.01 have 
been included in the calculation, this setting maybe too 
imprecise. Using data from [s^, different values for the 
coefficients have been obtained as c\ = 33.6, C2 = —249, 
C3 = 599, and C4 = —461. In this case, the critical /3 was 



FIG. 9: The full momentum dependence of the zero-mode 
electric (top panel) and magnetic (bottom panel) gluon prop- 
agator in four dimensions. To extend the temperature range 
outside the [—0.1,0.1] interval additional lattice sizes have 
been used, which are chosen to have an approximately con- 
stant number of lattice points as the 6 x 48^ lattices in this 
reduced temperature range. Dips along the t = line are 
artifacts from the graphical interpolation routine. 



determined more precisely as /3 = 2.2986, instead oi ^ — 
2.299, as used in the main part of the text. However, this 
leads only to a shift in t of the order of 0.01, i. e., reduced 
temperatures in the main text should be taken to have 
a systematic uncertainty of ±0.01. This accommodates 
the shift of the peak of the susceptibilities in figure S] 
compared to f = 0. Similar comments apply, of course, 
to the three-dimensional case. 



12 




FIG. 10; The electric screening mass and susceptibility deter- 
mined using a temperature-dependent renormalization pre- 
scription, as in the main text, and a temperature- independent 
renormalization prescription. 



2. Renormalization 

In four dimensions, the propagators have to be renor- 
mahzed. As a consequence, also the screening mass is 
rcnormalized by the same wave-function renormahzation 
factor. In the main text, this is done at each temperature 
individually at 2 GeV. Of course, since finite temperature 
is not introducing any new divergences [G^I , it is also pos- 
sible to renormalize once and for all at zero temperature. 
However, since the renormalization constants are both 
potentially volume-dependent and not available for all 13 
values investigated, this has not been done here. 

To assess the relevance of renormalization, it is shown 
in figure [TU] the effect of renormalizing only once at zero 



FIG. 11: The electric screening mass and susceptibility as a 
function of the reduced temperature at fixed spatial volume 
for two different volumes. 

temperature for the 4 x 46"^ lattice. Here, the /? depen- 
dence of the renormalization constants is ignored, given 
that it is expected to be only logarithmic in the scaling 
regime. The result show indeed no change within statis- 
tical errors. This was expected, as at a momentum of 2 
GeV temperature effects are still rather small. This may 
change at lager temperatures, but is not relevant for the 
present purpose. 



3. Finite volume effects 

In the main text, the spatial volume was not kept fixed, 
but small variations have been permitted to obtain a very 
dense mesh of temperatures around the critical one. In 
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FIG. 12: The electric screening mass as a function of the 
aspect ratio at T = 0.95Tc. The aspect ratios from left to 
right are, respectively, 1:4, 1:5, 1:5, 1:7, 1:7, and 1:9. 



particular, the spatial volumes at very low temperatures 
and very high temperatures are quite different. To assess 
the influence of this effect, the electric screening mass and 
susceptibility it shown at fixed spatial volume in figure 

m 

As is visible, the overall generic behavior over the 
whole temperature range is not significantly altered when 
using a fixed instead of a varying spatial volume. How- 
ever, as has already been observed in other calculations 
[2l| , it is visible that a finer lattice seems to increase the 
plateau directly before the phase transition. Of course, 
in the present case neither the temperature mesh, nor 
the spatial volumes, are exceptionally large. Therefore, 
this should be taken only as an indication, and requires 
further study. 



4. Aspect ratio 

As already indicative in the discussion of the finite 
volume effects, also the aspect ratio and the temporal 
discretization play a significant role. In figure [T^] the 
electric screening mass at a fixed temperature is shown 
for a number of different spatial volumes and aspect ra- 
tios. In the figure, the aspect ratio increases from left 
to right, approaching the infinite- volume and continuum 
limit such that the temperature extension remains finite, 
while the spatial extension diverges. It is visible that 
at small aspect ratios and volumes a significant depen- 
dence on these two parameters exist. However, towards 
the infinite- volume and continuum limit, this dependency 
quickly diminishes, if the limit is taken in the presented 
order, though the statistical error bars are rather large. 
The final value is in agreement with the one obtained 
in the main text. In contrast, at fixed aspect ratio and 



spatial volume a much larger imp act of going to the con- 
tinuum limit has been found jl9l . [2l| , indicating that the 
different approaches towards the desired limit have to be 
investigated carefully. 



Appendix C: Statistical errors 

The assessment of the quality of fits presented in the 
main text is a rather complicated problem. Usual 
checks can only be expected to give a reliable estimate if 
the fit form is linear, but the present one is non-linear, 
non-analytic, and non-continuously differentiable, and in 
the case of SU(3) not even continuous. To estimate the 
statistical error of the fit thus instead a comparison to 
a null hypothesis is used. This null hypothesis is the 
absence of any signal of the phase transition, i. e. , a trivial 
linear behavior 



Mf{t) =mQ + a{l + t). 



(CI) 



For a quantity not infiuenced strongly by the phase tran- 
sition this is an adequate description, as is shown in case 
of the magnetic screening mass in figure 1131 

Magnetic screening mass^ 

0.7 r 




FIG. 13: The magnetic screening mass on the 4 x 46 lattice 
close to the phase transition temperature compared to a linear 
fit. 

To estimate the error, the results for the electric screen- 
ing mass have also been fitted to the formula (jCip . Then, 
instead of using the average values for the electric screen- 
ing mass, a value within the 67% confidence interval for 
each temperature point was chosen such that it was clos- 
est to the fitted linear behavior. Then these points were 
fitted again within the fit ansatz ([T^. Since the ansatz 
([T2|) contains also this null hypothesis, this gives expo- 
nents as close as possible to the null hypothesis. To es- 
timate an opposite error, the fiuctuation were performed 
towards as far away as possible from the null hypothesis 
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(jCl|) . These two fits gave two further fit parameters. In 
the tables |T] and HI] the deviation of the fit parameters 
towards the null hypothesis is set as a subscript, while 
those as far away from the null hypothesis as possible 
are set as superscripts. These results have been used 
in the figures [3] and [S] In the four-dimensional case in 
the low-temperature phase the uncertainties in the fits is 
quite significant, a consequence of the compatibility with 



a constant within errors. Thus fits with free pre-factor 
and exponent can turn out to be optimal by combining 
very extreme combinations of both. The result should 
thus be rather interpreted as that the low-temperature 
behavior is within errors compatible with a constant. 

In the SU(3) case in figure [7] the error for the fit was 
determined in the same way. 
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